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Abstract 

We study the motion of an overdamped colloidal particle in a time-dependent non-harmonic 
potential. We demonstrate the first law-like balance between applied work, exchanged heat, and 
internal energy on the level of a single trajectory. The observed distribution of applied work is 
distinctly non-Gaussian in good agreement with numerical calculations. Both the Jarzynski relation 
and a detailed fluctuation theorem are verified with good accuracy. 

PACS numbers: 05.40.-a, 05.70.-a 



1 



Since more than a century, the first law relating the work applied to a system with 
both the exchanged heat and an increase in internal energy is one of the corner-stones 
of macroscopic physics. Its consistent formulation for a mesoscopic system like a driven 
colloidal particle, however, was suggested only about a decade ago (l| . Since on these scales 
thermal fluctuations are relevant, probability distributions for work, heat and internal energy 
replace the sharp values of their macroscopic counterparts. Various theoretical relations like 
the fluctuation theorem 0, jj, the Jarzynski relation [4], and the Hatano-Sasa relation ^] 
involving these distributions in different settings extend the second law to the mesoscopic 
realm at least as long as the notion of a constant temperature of the ambient heat bath 



remains meaningful (for a review, see [6]). Such theorems have been tested experimentally 
using both biomolecules manipulated mechanically [?} Q as well as colloidal particles in time- 



dependent laser traps [9 
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Common to all colloidal experiments, so far, is that these 
laser traps generate a harmonic potential albeit with a time- dependent center or "spring 
constant". Consequently, often the interesting distributions are Gaussian even though for 
certain quantities non-Gaussian distributions can occur 

In this Letter, we study the thermodynamics of single colloidal trajectories in a time- 
dependent non-harmonic potential which, generically, gives rise to non-Gaussian distribu- 
tions. Only for very short or very long trajectories, one expects Gaussian distributions even 
in this non-harmonic case |l3[. In particular, we identify applied work, exchanged heat 
and change in internal energy along a single trajectory and thus test the consistency of 
these notions on this level, or, put differently, illustrate the validity of the first law. We 
measure the distribution of work in the non-Gaussian regime and compare it to theoretical 
prediction. Such a comparison does not involve a single fit parameter since all quantities are 
measured independently, which is another advantage of colloidal systems. Finally, we test 
the Jarzynski relation which expresses the free energy difference between two equilibrium 
states in terms of the nonequilibrium work spent in the transition between the two states. 
Such an illustration of the Jarzynski relation in the non-Gaussian regime comes timely given 



ongoing theoretical criticism of its validity 
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In our study, particle trajectories were determined using total internal reflection mi- 
croscopy (TIRM), where a single colloidal particle is illuminated under evanescent field 
conditions. This field is created by total internal reflection of a laser beam at a glass-water 
interface. The scattered intensity of a bead near the interface is proportional to exp(— £ z), 
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with the decay length of the evanescent field and z the particle-wall distance Q|. Mea- 
suring the scattered intensity of a fluctuating Brownian particle as a function of time thus 
yields its vertical position with a spatial resolution of about 5 nm. 

We used an aqueous suspension of highly charged polystyrene beads with radius R = 
2 /im, which were illuminated with light of wavelength A = 658.5 nm. The particle concen- 
tration was sufficiently low to guarantee that there was only a single particle within the field 
of view. The penetration depth was adjusted to ~ 200 nm and the scattered intensity 
was monitored with a photomultiplier at a data acquisition rate of v = 2 kHz. An additional 
focused laser beam (A = 1064 nm, power P ~ 2 mW) was directed vertically from the top, 
which confined the particle motion to an one- dimensional trajectory in ^-direction. 

To drive the colloidal particle between two equilibrium states, it was subjected to the light 
pressure of another optical tweezers (A = 532 nm, beam waist about 17 /im, P < 60 mW), 
which was incident from below into the sample cell (lower tweezers, see inset of Fig. ^) . The 
intensity of this laser beam was varied with an electro-optical modulator connected to a 
computer-controlled waveform generator. We modulated the laser intensity according to a 
time-dependent symmetric protocol I(t) = I(t s — r) in the time interval < r < t s , where 
t s is the pulse duration (see Fig. |2J) |25[. To ensure that the system is out of equilibrium, 
t s must be smaller than the particle relaxation time t r ~ 480 ms. On the other hand, when 
repeating the experiment the pause t p between two consecutive pulses must be longer than 
t r to guarantee equilibration of the system. To meet both conditions in our experiments we 
have chosen t s = 120 ms and t p = 700 ms. 

The total time- dependent potential acting on the particle at distance z is given by [lj| 

V(z, t) = A exp(-nz) + B z + C I(t)z. (1) 

The first term describes the double-layer interaction between the negatively charged col- 
loidal particle and the likely charged wall with A depending on the corresponding surface 
charges and k _1 the Debye screening length, which depends on the salt concentration in the 
suspension. The second term accounts for the weight of the particle and the additionally 
exerted light pressure from the upper tweezers, which both depend linearly on the particle 
distance z j^. The last term considers the time-dependent optical forces induced by 
the lower tweezers. Experimentally, particle-wall potentials are easily obtained by measur- 
ing the distance probability distribution p(z) of a colloid in front of a wall. In thermal 
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FIG. 1: (color online) Particle wall interaction potentials for three different intensities of the lower 
optical tweezers [decreasing power from (1) to (3)]. The solid line shows the fit according to Eq. (fTJ) 
with k^ 1 ~ 25 nm. Inset: Light pressure versus tweezers intensity. The light pressure is a linear 
function of the laser intensity. 

equilibrium, i.e. for J(r) = const., the potentials are given by V(z) = —k B T\np(z) up 
to a constant. Here, T is the temperature of the environment and k^T is the Boltzmann 
constant. The symbols in Fig. Q show V(z) obtained for three different intensities of the 
lower tweezers. The solid line is a fit to Eq. ^ (exemplarily shown only for one data set) and 
clearly demonstrates that the particle is moving in a non harmonic potential. It can be seen 
in Fig. ^ that the light pressure of the lower tweezers reduces the slope of the linear part of 
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171 ] . The inset of Fig. [T] shows 



the potential, as already demonstrated by other authors 
the expected linear dependence on the intensity of the light pressure. 

Fig. 121 shows exemplarily the trajectory of a particle driven by the time dependent po- 
tential While during the first pulse the particle is strongly forced towards the surface, 
thermal fluctuations support it to move against the applied force away from the wall during 
the second pulse. This clearly demonstrates that the particle is strongly coupled to the 
surrounding heat bath. 

In _ to recent expenntenta, stndies 333130, we want to test elementally 
to which precision energy conservation of the particle on its trajectory is maintained. This 
does not only provide a rigorous check of the experimental technique, data analysis, and the 
energy resolution but also demonstrates the interplay of applied work and exchanged heat 
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FIG. 2: (color online) Measured tweezers intensity and particle trajectory. During the first pulse 
the particle is pressed towards the surface. During the second pulse thermal fluctuations support 
the particle and it is able to move away from the wall. Hence the applied work is positive for the 
first pulse and negative for the second. 

when the system is non-adiabatically driven. Therefore, in addition to the work W exerted 
on the particle we need to determine its heat exchange Q with the environment. 
The Brownian motion of the colloidal bead is described by the Langevin equation 
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dV , 

oz 



(2) 



with 7 the friction coefficient and £ the stochastic force. According to jlj the incremental 
change of heat dQ and work dW is then given by 

dV 



dV 

dW = -r—dr, dQ 
or 



dz 



■dz. 



(3) 



Integration along a single trajectory z(t) then leads to the work functional 
W[z{r)\ - 







Or 



Ti Z(Tj, 



(4) 



where we have used Eq. (JTJ). The right hand side of Eq. (j^J) accounts for the discrete sampling 
of the particle trajectory during our experiments with rate v — i at times Ti = idt. Because 
the velocity autocorrelation of a Brownian particle decays on a timescale of some 10 ns, the 
velocity 

rn+i 

z(n) —v\ dr z(r) (5) 
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FIG. 3: (color online) (a) The quantities — Q, W and AV for about 100 periods of the protocol 
J(r). (b) Distribution histogram of 5 = W — Q — AV, the experimentally observed "deviation" 
from the first law of thermodynamics. 

determined from a trajectory measured with v = 2 kHz is not identical to the instant particle 
velocity z. However, since dV/dz varies on a timescale much larger than 5t, the heat along 
a single trajectory z(r) can be written as 



With the above sign convention the heat Q is negative (positive) when extracted (delivered) 
from (to) the thermal environment. Since we have full knowledge of the time-dependence of 
I, V, and z, both quantities W and Q can be determined from our experiments. 
Introducing AV = V(z(t s ),t B ) — V"(^(0),0) we finally obtain a stochastic version 



of the first law of thermodynamics. Fig. a) shows work W, heat Q and change of inner 
energy AV for the trajectory of a single particle where the protocol 7(r) was repeated about 
100 times. For W and Q maximal energies of about 15 k^ T are observed, whereas AV is 
on the order of a few T. Obviously, Q and W are not independent quantities. Usually 
trajectories resulting in a large work W are also accompanied by a large value of Q. But 
only when taking all three energies in Eq. (JIJ) into consideration, the distribution of the 
deviation shown in Fig. Efb) is centered around zero, having a half-width of about 0.7 k^T. 
Assuming that the three terms have the same contribution to the total error, the energy 
error of these energies is about one quarter of k B T. 
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W[z(r)] -Q[z(t)] -AV = 



(7) 
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FIG. 4: (color online) Non-Gaussian work distribution. The data was taken from about 16000 
trajectories, where the average work done on the particle was about 2.4 T. The solid line shows 
the Pearson type III distribution 18] corresponding to the theoretically obtained moments. Inset: 
Logarithm of the ratio of the probability to find trajectories with work — W to those with work 
+W . The solid line shows the expected curve ©• The deviation is due to the poor statistics of 
large negative work values W < — 4 fee T. 

The measured work distribution in Fig. H] is distinctly non-Gaussian and exhibits an 
asymmetry towards larger work values. Whereas the first law is valid along a single trajectory 
as demonstrated, fluctuation theorems considering probability distributions can be regarded 
as an extension of the second law. This becomes evident when looking at the Jarzynski 
relation ^| 



-W/ks T 



> 



-AF/k B T 



(8) 



where AF is the change of free energy between two equilibrium states and the brackets 
(■ • •) represent the average over the work distribution spent in a transition between these 
states. Eq. (JSJ) immediately leads to (W) > AF, a formulation of the second law for driven 
systems on the mesoscopic scale. A test of the Jarzynski relation using the data shown in 
Fig. EJyields (exp(— W/k^ T)) ~ 1.03 in agreement with AF = for the symmetric protocol 
I(t). In addition to the integral theorem (jHJ) we also test the somewhat stronger detailed 
fluctuation theorem 



P{-W)/P{+W) = e 



-W/ka T 



(9) 



which holds for time-symmetric protocols J(r) = I(t s — r) (see inset of Fig. HJ |l9l l20l|. Here 
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mi [feT] 


m 2 [(^T) 2 ] 


m 3 [(^T) 3 ] 


Exp. 


2.4 


11.7 


67.8 


Theo. 


2.4 


11.6 


63.7 



TABLE I: Comparison between theoretically predicted and measured moments of the work prob- 
ability based on the data shown in Fig. 01 

the probability P(— W) that a negative work value occurs is compared to the probability 
P(+W) of a positive value of same magnitude. 

In order to compare the measured histogram in Fig. 0] to the theoretical prediction we 
calculate the probability distribution solving the Fokker-Planck equation |l3t \21\ 



dp dV dp 

dt dr dw 



(10) 



Here, p(z, w, t) is the joint probability of the particle to be at time t a distance z away from 
the wall and to have accumulated an amount of work w up to this time. The Fokker-Planck 
operator j^] 

c _3 (D ± dV d\ 

governs the dynamics of the particle, where 

D L {z)^D,[l + R/{z-R)}- 1 (12) 

is the diffusion coefficient for perpendicular motion near a surface j^J. The free diffusion 
constant for a particle with radius R = 2 pm is Dq ~ 0.1/xm 2 /s at room temperature 
T = 293 K. Since we start in equilibrium with no work spent on the particle yet, the initial 
distribution needed to solve Eq. (fTTH) is 

exp(-V(z,0)/k B T) 



p(z, w, 0) = S(w)- 



(13) 



/ dz exp(-V(z,0)/k B TY 

Eq. (flT)|) is a Fokker-Planck equation in two space dimensions z and w including a singular 
initial condition, which numerically is not easy to handle. We therefore multiply Eq. (|10J) 
with w n and integrate over w. After one integration by parts we obtain the inhomogeneous 
evolution equation 

(14) 



dt 



dV 

CM n + n— M n _! 
dr 



S 



for the conditional nth moment of the work 

/+oo 
dw w n p{z,w,t). (15) 
-oo 

The actual nth moment m n then follows simply by integrating over z and thus adding 
the contributions of all possible final positions of the particle. The function M (z,t) is the 
probability distribution of the position z of the particle and hence mo = 1. Table H] compares 
the numerically and experimentally obtained first three moments of the data shown in Fig. 0] 
We stress that this good agreement does not involve a single fit parameter. 

In summary, we have confirmed experimentally both a stochastic formulation of the first 
law and various recent theoretical ramifications of the second law in a time-dependent non- 
harmonic potential, where the underlying distributions are typically non-Gaussian. In the 
next step, non-harmonic systems with broken detailed balance should be investigated to test 

nn 

theorems |5j, |24| which, so far, have been under experimental scrutiny in the harmonic case 
only Q. 
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